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O . Abstract 

CN ' Supersymmetry (SUSY) with a long-lived stau is an attractive scenario in the LHC 

\Q • experiments because one can directly observe stau tracks in each SUSY event, and 

, thus precise measurements of SUSY particle masses are possible. In this scenario, 

we discuss the possibility to observe/measure parity violation in interactions among 
SUSY particles. Such a measurement will be important in determining spins and 
chiralities of SUSY particles. We use the last step of the cascade-decay chain: — )• 
rS \ TT —7- f(lvv\ where the polarization of the tau lepton can be determined statistically 

■ by looking at the energy distribution of the final state lepton. Comparing with the 

theoretical formula of the neutralino differential decay width, one can extract the size 
of parity violation in the interaction vertices among the stau, the tau lepton and the 
neutralino. We perform a Monte Carlo simulation to see if the effect is visible at the 
LHC experiments. 
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1 Introduction 



If supersymmetry (SUSY) is the solution to the hierarchy problem, the superpartners should 
weigh around a few hundred GeV to TeV energy range, which is accessible at the LHC 
experiments. The way to discover SUSY or to measure the properties of SUSY particles are 
quite different depending on the pattern of the superparticle spectrum. In particular, the 
property of the lightest SUSY particle (LSP) is important for the search strategies. 

Motivated by an explanation of dark matter of the Universe, many collider studies of 
SUSY models have assumed the case with the LSP being the neutralino (admixture of the 
Higgsinos and the gauginos). In the neutralino LSP case, the final states of the SUSY events 
always contain invisible neutralinos, and one needs various measurements to be combined in 
order to determine sparticle masses. The situation drastically changes if one assumes that 
one of the charged sleptons is lighter than the lightest neutralino. Such a scenario is possible 
if the slepton is not absolutely stable. For example, the slepton can decay into a lepton and 
gravitino if kinematically allowed. In particular, the scalar tau lepton (stau) can easily be 
lighter than the neutralino in many SUSY models due to the quantum correction through 
the large Yukawa interaction. 

If the life-time of the stau is longer than the time scale of the collider experiments (a 
few nsec), we will see charged tracks left by staus in each SUSY events, instead of the 
missing momentum in the neutralino LSP scenario. Methods to search for a long-lived 
charged particle in hadron colliders have been studied in Refs. [H [21 El IH E] • Also, various 
advantages in studying SUSY models at the LHC have been reported. The stau momenta 
and velocities can be measured by analyzing the stau tracks, from which the stau mass can be 
extracted with a good accuracy [6l[71|8]. The measurements of other sparticle masses [9l [TOl 
[TT| [T2I IT3l [H] and the spin measurement of the stau [I5] have also been discussed. Methods 
to measure the stau hfe-time have been proposed in various contexts (TBI HJl [TSl UHl 120] . 
The chargino-neutralino production and subsequent decay processes at the LHC have been 
studied in Ref. and the possibilities to observe parity/CP violations have been discussed. 

In this paper, we discuss the possibility to measure the tau polarization in the neutralino 
decay, x° ~^ t't, where the neutralinos are mainly produced by the cascade decays of colored 
SUSY particles. The tau polarization carries information on whether the lightest stau is the 
partner of the right- or left-handed tau lepton, and also whether the decaying neutralino is 
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mainly composed of the gauginos or the Higgsinos. The polarization can be determined by 
looking at the energy distribution of the r-decay product [22] due to the fact that the weak 
interaction maximally violates parity. We use the leptonic decay mode of the tau lepton, 
r — !■ luu, in this work. We first show the formula for the lepton-energy distribution in the 
neutralino decay, and by using it we fit the data from the Monte Carlo simulation. We find 
that the size of parity violation (the tau polarization) can be measured in a simple model 
where only one of the neutralinos contributes to the final state. In a more complicated 
case, where there are two neutralinos contribute to the same final state, we can correctly 
reproduce the signs of parity violation in each decay vertex. 

2 Neutralino decay 

In this section, we present the differential decay width of ^ ^ ^{^i^i^r) as a function of 
the lepton-energy fraction. The parity asymmetry of this decay process carries information 
on the chirality of the stau and the composition of the neutralino, i.e., whether it is gaugino- 
like or Higgsino-like. 

The measurement of the tau polarization has also been studied in the neutralino LSP 
case. The effects of the tau polarization in the di-tau invariant mass distribution in cascade 
decays were studied in Refs. [231 1211 125] • Ref . |2S] studied the px distribution of the softest r- 
jets to extract the tau polarization in the co-annihilation region in the mSUGRA model. The 
tau polarization measurements at e^e~ colliders were studied in the stau-pair productions 
in Refs. [271 |2S1 122] • The measurements of CP violation at e"'"e~ colliders were also studied 
in Refs. [301 [3ll |32] . 

2.1 Decay width 

The polarization information of the tau lepton is imprinted to the energy distribution of the 
lepton in the leptonic tau decay. We show in this subsection the formula for the distribution 
of the lepton-energy fraction in the neutralino decays. 

The relevant interaction Lagrangian for the neutralino decay is 

£ = ^(^?l,Pl + 9mPR)rf^ + H.c. , (1) 
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where Pl and Pr are the chirahty projection operators, and ql and are couphng constants. 
We will discuss the relation to the underlying model parameters in the next subsection. The 
tau lepton is polarized if there is parity violation, i.e., \gL\ ^ \gR\- The exphcit calculation 
shows that the differential decay width of the process, 'X^i ^ tt ^ {lviVr)f., is 

L^= (1-^) (5 + 5^-4^2) _aW(i + ^_ 8^2) dz , (2) 

f lept -J L 

where Fiept is the partial decay width of the leptonic mode [22]. The variable z is the 
lepton-energy fraction: 

where the energy of the lepton and r [Ei and Et) are those in the rest frame of The 
formula ([2]) is independent of the tau charge or the flavor of the final state lepton, i.e.., 
electron or muon. The parameter af^ represents the tau polarization and it is expressed in 
terms of the Lagrangian parameters: 

ii) _ \9Li? - \9Ri? 

^ ~ \n |2 I 1^ 12 • W 

\gLi\ + \gRi\ 

The z distributions in Eq. ([2]) for = — 1,0,+1 are shown in the left panel of Fig. [H 
The value = — 1(+1) corresponds to the decay into right (left) handed r. The lepton 
tends to be emitted to the same direction of r if oat = 1 (left-handed r). 

2.2 Coupling constants and the model parameters 

The coupling constants in Eq. ([1]) are expressed in terms of the neutralino- and the stau- 
mixing parameters as follows: 

gii = {--^Oii - ^0^2) sin6'f ^^^Oi3Cos6'f , (5) 

V V2 V2 / V cos 13 



9Ri 



-OissinOf + \/2giOii cos Or . (6) 



V cos (3 



where gi is the coupling constant of the U(1)y gauge interaction, g2 is that of SU(2)l, tan/3 
is the ratio of the two vacuum expectation values (VEV) of the Higgs fields and v is the 
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dr/ Tieptdistribution dr/ Tieptdistribution with a model of lepton Pj cut 
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Figure 1: Distributions of the lepton-energy fraction {z = Ei/Er) in the neutrahno decay. 
The z dependence in Eq. with otv = 0, ±1 are shown in the left panel. The right panel 
is the same figure with a simple model to include the effects of the lepton px cut. 



VEV of the Higgs field: 



tan /3 

2 



V 



(7) 
(8) 



The matrix elements Oij are those of a unitary matrix which diagonalizes the neutralino 
mass matrix, 

/Ml 
M2 



-^cos/3 



giv 



sm 



-2ifcos/5 ^cos/3 



V2 

-ffsin/3 



y^sin/3 -^sin/3 



-/i 




such that 



OM^oO' 



The mixing angle of the stau is defined by 

Ti = tl sin 9f + cos 9f 



(9) 



(10) 
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For example, let us consider the case in which the lightest neutralino (x?) and the lightest 
stau (fi) are mainly the Bino and tr, respectively, i.e., cosOf ~ 1 and On ^ O12, O13. In 
this case, the largest contribution to the x? decay is the second term in Eq. ([6]). Therefore, 
the parity violation parameter a"^ in Eq. (|2]) is a^^ ~ — 1. 

In the case where the decaying neutralino is mainly the Wino and the lightest stau is tr, 
the situation is different. Since the pure Wino does not couple to the right-handed stau, the 
decay occurs either through the neutralino mixing or the stau mixing. In both cases, the 
parity violation is a^v > in contrast to the previous example. 

3 Monte Carlo simulation 

In this section, we demonstrate by using a Monte Carlo simulation that the distribution 
obtained in the previous section is observable at the LHC experiments. 

3.1 Basic set up 

We generated events of SUSY particle productions at a pp collider at ^/s = 14 TeV (LHC) 
by using the HERWIG package [33]. We use CTEQ5L [M] for the parton distribution function. 
The TAUOLA package [35] is used for the r decays and the events are passed through the 
AcerDET detector simulator [36] where the lepton momentum are smeared to simulate the 
detector effects. 

The final state contains two stau tracks. The identification of those staus can be used to 
eliminate the background from the Standard Model processes. In order to distinguish staus 
from muons, we impose following selection cuts on the candidate stau tracks: 

• 0.40 < (3 < 0.91 , 

• Pt > 10 GeV , 

• < 2.4 . 

Here and hereafter, we do not take into account the momentum and the velocity resolutions 
of the staus, i.e., the parton level information is used. We assume in the following that the 
stau and neutralino masses are known. It has been reported in Refs. [3 El [9|, [TOl [TTl Hlj that 
those quantities are measurable with a good accuracy. 
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3.2 Lepton px cut and the deformation of the z distribution 

We also require pt > 15 GeV for the lepton momentum. This cut affects the shape of the z 
distribution in Eq. ([2]). We model the effect of the pt cut in the z distribution by multiplying 
the following weight: 

^Erfcf--^) , (12) 



2 V V2a^ 



where p™" is 15 GeV, and C and parameters to be determined by fitting the data. 

We plot in the right panel of Fig. [1] the z distribution after multiplying this factor with 



C = 150 GeV and = 6.0 GeV. 



3.3 Model I 



We use the MSSM without flavor mixing and ignore the Yukawa interactions of the first and 
second generations. As the first example, we take a simple model in which only one of the 
neutralinos and one of the staus are significantly lighter than others, and the lightest 
neutralino and the lightest stau are almost the Bino and right-handed, respectively. The 
value of the parity violation parameter is 

= -0.99 , (13) 

in this case. 

The MSSM parameters we take are listed in Tab. [H We follow the convention of SUSY 
Les Houches Accord [37j. The x° and f masses are calculated to be 

m^o = 195 GeV , m^ = 118 GeV . (14) 

We used the ISA JET package [38] for the calculation of the mass spectrum and the branching 
ratios. We generate 10,000 SUSY events with the parameter set given in Tab. [H The number 
of events corresponds to the integrated luminosity of 14.5 fb~^. 

The neutralinos are mainly produced by the cascade decays of the squarks and gluinos. 
The value of z (the lepton-energy fraction) in each candidate event can be calculated from 
the measured invariant mass of the lepton and the stau. 



(15) 
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Parameter 


Value 


Ml 


195.7 [GeV] 


M2 


800.0 [GeV] 


M3 


1003.7 [GeV] 


tan P 


39.46 


u 


800.0 [GeV] 


B 


694 [GcV] 


Mo, , Moo 


1400 [GeV] 


Ml, , Ml, 

i^l J -'-'2 


518[GeV] 

L J 


Mq, 


1241 [GeV] 


Mls 


487[GeV] 


Muc 


1321 [GeV] 




13 11 [GeV] 


Mcc 


1321 [GeV] 


Msc 


1311 [GeV] 


Mec 


278 [GeV] 


M^c 


277[GeV] 


Mfec 


1210 [GeV] 


Mjc 


1075 [GeV] 


M^c 


160[GeV] 




-503 [GeV] 


A;. 


-567[GeV] 


A. 


-48.3[GeV] 



Table 1: Parameters of Model I. 
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Figure 2: The histograms of the t-1 invariant mass (left) and z (= Ei/E^-) in Model I. 



where Mj^i is 

M',i = {p,+pi)\ (16) 



We show in Fig. [2] the histograms of the invariant mass (left) and the z parameter 
(right). In these plots, we select all possible combinations of a stau and a lepton with the 
opposite charges in each candidate event. In order to eliminate the contributions from wrong 
combinations and from background events such as leptons from IV-boson decays, we subtract 
the same distributions calculated by using the combinations with the same signs [H]. After 
the selection cuts and the charge subtraction explained above, 978 events remained. 

By fitting the z distribution in Fig. [2] with the function in Eq. ([2]) multiplied by a factor 
in Eq. f[T^ . we obtain the ajy parameter to be: 

= -0.95 ±0.32 . (17) 

Comparing with the input value in Eq. (1131) . we can see that the parity violation is measured 
quite successfully. 
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Figure 3: The histogram of the t-1 invariant mass in Model II. 



3.4 Model II 

As the second example, we take a model where the Wino-like neutralino is not so heavy 
compared to the Bino-like one as is so in the models motivated by the GUT relation of the 
gaugino masses. The model parameters are listed in Tab. [2J The neutralino and stau masses 
are 

m^o = 195 GeV , m^o = 359 GeV , = 118 GeV . (18) 
The parity asymmetries a^'' and calculated to be 

= _o.99 , = 1.0 . (19) 

Unlike the case of Model I, we cannot reconstruct z in the event-by-event basis since we 
do not know whether the decaying neutralino is Xi oi' X2- Therefore, we need to directly fit 
the invariant mass distribution which contains events of both neutralinos. 

We generated 90,000 SUSY events with this parameter set. The number of events corre- 
sponds to the integrated luminosity of 104.6 fb~^. We show the invariant mass distributions 
(Fig. [3]) using the same technique as we used for the Model I. After the cuts and the sub- 
traction explained before, 7,677 events remained. We can see contributions from decays of 
two kinds of neutralinos. 
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Parameter 


T r 1 

Value 


Ml 


195.7 [GeVJ 




363.9 [GeVJ 


Ms 


1003.7 [GeV] 


tanp 


39.46 


1^ 


800.0 [GeVJ 


B 


694 [GeVJ 


71 /T 71 /T 

Mq„Mq^ 


1 A f\f\ \ T rl 

1400 [GeVJ 


Mr Mt 


518fGeVl 




1241 [GeVJ 




487[GeVj 




1321 [GeVJ 




13 11 [GeVJ 




1321 [GeVJ 




1311 [GeVJ 




278[GeVj 




277[GeVj 




1210[GeVj 


Mfc 


1075 [GeVJ 


Mrc 


160 [GeVJ 


At 


-504 [GeV] 


At, 


-567[GeVj 


Ar 


-48.3[GeVj 



Table 2: Parameters of Model II. 
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Before fitting the distribution in Fig. [31 we cliecked the validity of Eq. ([2]) combined 
with Eq. f|T2|) for the Model II. We did it by separating the Xi ^"^^ the X2 events by using 
information from the event generator and performed the fitting for each event set. Reasonable 



values for and obtained: 



(1) .^A .(2) 

(1) _ _n sn -I- n 1 Q 



a)^^ = -0.80 ±0.13 , a^^ = 0.94 ±0.084 . (20) 

For the fitting of the histogram in Fig. [31 we use the sum of two fitting functions corre- 
sponding to each neutralino: 

^"P I' ^ \ 1 / ^min r~< ^ \ 

(21) 

where p™" = 15 GeV, and dFj/rjgp^ is the function in Eq. ([2]). The energy fraction Zi is 
related to the f-l invariant mass, Mfi, as 




2 2 



(22) 



There are eight fitting parameters a!~^ , Ni, Ci and o"*^ for z = 1, 2 in the function. 

Since it is difficult to find the best fit for such many parameters, we here make a (semi) 
theoretic assumption on a relationship between two parameters Ci and C2 in order to reduce 
the number of parameters. In Eq. f[2T]) . the combinations of CiZi represent the average values 
of from the Xi decay in the laboratory frame. Since the neutralinos are produced mainly 
from heavy colored particles, it is likely to be highly boosted and thus the tau lepton from 
the neutralino decay is pointing to the similar direction to the one of the neutralino in the 
laboratory frame. Therefore, there is an approximate relation: 



(P^) = ^^{PT) ^ ( 1^ ) ^) ' (23) 







where {Er/E^o) is 



,2 _ _2 

Xi 



^ ^ ' (24) 



E / 2m\ 
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The average transverse momentum of the neutrahno, {p^)^ is independent of the neutrahno 
mass if it is highly boosted. Therefore, the ratio of C1/C2 is estimated to be 



^0 m^o 

X2 Xi 



C2 mK m\ - ml 

Xl X2 



(25) 



By imposing the above relation, we could fit with the function in Eq. (|2T]) : 

= _o.66±0.27 , aSJ^ = 0.63 ±0.13 . (26) 

We obtain the correct signs for each parity violation although the magnitude of the par- 
ity violation is obtained to be smaller than the theoretical input values (—0.99 and 1.0, 
respectively) due to our simple ansatz on the fitting function. 



4 Summary 

We studied the polarization measurement of the tau lepton at the LHC in the long-lived stau 
scenario. The polarization of the tau lepton from the neutrahno decays carries information 
on whether the neutrahno is Higgsino-like or gaugino-like and the chirality of the partner of 
the stau. 

We have shown that the r polarization can be measured by fitting the distribution of 
the lepton-energy fraction [Ei/Et-) in the leptonic r decays. We performed Monte Carlo 
simulations for two parameter sets of the MSSM. The first example we took is a simple 
case where the Bino-like neutrahno is significantly lighter than others. In this case, we can 
directly fit the energy fraction by the theoretical function and the parity violation can be 
measured successfully. 

In the second case where the Wino-like neutrahno also contributes to the same final state, 
we fit the f-l invariant mass distribution with the two contributions summed. Although 
the fitting gives a milder asymmetry compared to the theoretical inputs, we can obtain 
qualitatively correct features of parity violation in each neutrahno decay. 
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